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Abstract 

We study transmission stability and dynamics of pulse amplitudes in 7V-channel 
soliton-based optical waveguide systems, taking into account second-order dis¬ 
persion, Kerr nonlinearity, delayed Raman response, and frequency dependent 
linear gain-loss. We carry out numerical simulations with systems of N coupled 
nonlinear Schrodinger (NLS) equations and compare the results with the pre¬ 
dictions of a simplified predator-prey model for Raman-induced amplitude dy¬ 
namics. Coupled-NLS simulations for single-hber transmission with 2 < A < 4 
frequency channels show stable oscillatory dynamics of soliton amplitudes at 
short-to-intermediate distances, in excellent agreement with the predator-prey 
model’s predictions. However, at larger distances, we observe transmission 
destabilization due to resonant formation of radiative sidebands, which is caused 
by Kerr nonlinearity. The presence of linear gain-loss in a single hber leads to 
a limited increase in transmission stability. Signihcantly stronger enhancement 
of transmission stability is achieved in a nonlinear A-waveguide coupler due to 
efficient suppression of radiative sideband generation by the linear gain-loss. As 
a result, the distances along which stable Raman-induced dynamics of soliton 
amplitudes is observed are signihcantly larger in the waveguide coupler system 
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compared with the single-fiber system. 
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1. Introduction 

Transmission of information in broadband optical waveguide links can be 
significantly enhanced by launching many pulse sequences through the same 
waveguide mill [HIS]. Each pulse sequence propagating through the waveguide 
is characterized by the central frequency of its pulses, and is therefore called 
a frequency channel. Applications of these multichannel systems, which are 
also known as wavelength-division-multiplexed (WDM) systems, include fiber 
optics transmission lines muins], data transfer between computer processors 
through silicon waveguides miiiiHi, and multiwavelength lasers m nni HU m]. 
Since pulses from different frequency channels propagate with different group 
velocities, interchannel pulse collisions are very frequent, and can therefore lead 
to error generation and cause severe transmission degradation mm mu 13 [13 

E]. 

In the current paper, we study pulse propagation in broadband multichannel 
optical fiber systems with N frequency channels, considering optical solitons 
as an example for the pulses. The two main processes affecting interchannel 
soliton collisions in these systems are due to the fiber’s instantaneous nonlinear 
response (Kerr nonlinearity) and delayed Raman response. The only effects of 
Kerr nonlinearity on a single interchannel collision between two isolated solitons 
in a long optical fiber are a phase shift and a position shift, which scale as l/A/3 
and l/A/3^, respectively, where A/3 is the difference between the frequencies of 
the colliding solitons USUI US]. Thus, in this long fiber setup, the amplitude, 
frequency, and shape of the solitons do not change due to the collision. However, 
the situation changes, once the finite length of the fiber and the finite separation 
between the solitons are taken into account m In this case, the collision 
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leads to emission of small amplitude waves (continuous radiation) with peak 
power that is also inversely proportional to A/?. The emission of continuous 
radiation in many collisions in an A-channel transmission system can eventually 
lead to pulse-shape distortion and as a result, to transmission destabilization 
m- The main effect of delayed Raman response on single-soliton propagation 
in an optical fiber is an 0{efj) frequency downshift, where cr is the Raman 
coefficient [Min [2D]. This Raman-induced self frequency shift is a result of 
energy transfer from high frequency components of the pulse to lower frequency 
components. The main effect of delayed Raman response on an interchannel 
two-soliton collision is an 0(e_R) amplitude shift, which is called Raman-induced 
crosstalk It is a result of energy transfer from 

the high frequency pulse to the low frequency one. The amplitude shift is 
accompanied by an 0(ej{/A/3) collision-induced frequency downshift (Raman 
cross frequency shift) and by emission of continuous radiation [5S1 US] HH 
nziiiniiiEi. Note that the Raman-induced amplitude shift in a single collision is 
independent of the magnitude of the frequency difference between the colliding 
solitons. Consequently, the cumulative amplitude shift experienced by a given 
pulse in an A-channel transmission line is proportional to a result that is 
valid for linear transmission imiiiiisniiDii, conventional soliton transmission 
[SSIEllESlEZ], and dispersion-managed soliton transmission [26]. Thus, in a 
100-channel system, for example, Raman crosstalk effects are larger by a factor 
of 2.5 X 10^ compared with a two-channel system operating at the same bit rate 
per channel. For this reason, Raman-induced crosstalk is considered to be one of 
the most important processes affecting the dynamics of optical pulse amplitudes 

in broadband fiber optics transmission lines [I1I2I1I11I21ISI1I31I31IM1I35]. 

The first studies of Raman crosstalk in multichannel hber optics transmis¬ 
sion focused on the dependence of the energy shifts on the total number of 
channels [21| . as well as on the impact of energy depletion and group velocity 
dispersion on amplitude dynamics isniiM]. Later studies turned their attention 
to the interplay between bit-pattern randomness and Raman crosstalk in on-off- 
keyed (OOK) transmission, and showed that this interplay leads to lognormal 


3 


statistics of pulse amplitudes m El HZl in EH ISZ]- This finding means that 
the nth normalized moments of the probability density function (PDF) of pulse 
amplitudes grow exponentially with both propagation distance and n^. Further¬ 
more, in studies of soliton-based multichannel transmission, it was found that 
the nth normalized moments of the PDFs of the Raman self and cross frequency 
shifts also grow exponentially with propagation distance and n^ |33l [34] . The 
exponential growth of the normalized moments of pulse parameter PDFs can 
be interpreted as intermittent dynamics, in the sense that the statistics of the 
amplitude and frequency is very sensitive to bit-pattern randomness [33| |34l |38] . 
Moreover, it was shown in Refs. |53J|331|3S] that this intermittent dynamics has 
important practical consequences in massive multichannel transmission, by lead¬ 
ing to relatively high bit-error-rate values at intermediate and large propagation 
distances. Additionally, the different scalings and statistics of Raman-induced 
and Kerr-induced effects lead to loss of scalability in these systems |55] . 

One of the ways to overcome the detrimental effects of Raman crosstalk 
on massive OOK multichannel transmission is by employing encoding schemes, 
which are less susceptible to these effects. The phase shift keying (PSK) scheme, 
in which the information is encoded in the phase difference between adjacent 
pulses, is among the most promising encoding methods, and has thus become 
the focus of intensive research [39l|4D]. Since in PSK transmission the infor¬ 
mation is encoded in the phase, the amplitude patterns are deterministic, and 
as a result, the Raman-induced amplitude dynamics is also approximately de¬ 
terministic. A key question about this deterministic dynamics concerns the 
possibility to achieve stable steady-state transmission with nonzero predeter¬ 
mined amplitude values in all channels. In Ref. |30j . it was demonstrated that 
this is not possible in unamplified optical fiber lines. However, the experiments 
in Refs. [miig showed that the situation is very different in amplified multi¬ 
channel transmission. More specifically, it was found that the introduction of 
amplification enables transmission stabilization and significant reduction of the 
cumulative Raman crosstalk effects. In Ref. |43| . we provided a dynamical ex¬ 
planation for the stabilization of PSK soliton-based multichannel transmission. 
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by demonstrating that the Raman-induced amplitude shifts can be balanced 
by an appropriate choice of ampliher gain in different channels. Our approach 
was based on showing that the collision-induced dynamics of soliton amplitudes 
in an A^-channel system can be described by a relatively simple A'^-dimensional 
predator-prey model. Furthermore, we obtained the Lyapunov function for the 
predator-prey model and used it to show that stable transmission with nonzero 
amplitudes in all channels can be realized by overamplification of high frequency 
channels and underamplification of low frequency channels. 

All the results in Ref. [43] were obtained with the A^-dimensional predator- 
prey model, which is based on several simplifying assumptions, whose validity 
might break down with increasing number of channels or at large propagation 
distances. In particular, the predator-prey model neglects high-order effects due 
to radiation emission, intrasequence interaction, and temporal inhomogeneities. 
These effects can lead to pulse shape distortion and eventually to transmis¬ 
sion destabilization (see, for example, Ref. [H])- The distortion of the solitons 
shapes can also lead to the breakdown of the predator-prey model description 
at large distances. For example, the relation between the onset of pulse pat¬ 
tern distortion and the breakdown of the simplihed model for dynamics of pulse 
amplitude was noted (but not quantified) in studies of crosstalk induced by 
nonlinear gain or loss [44l HU |46] . In contrast, the complete propagation model, 
which consists of a system of N perturbed coupled nonlinear Schrodinger (NTS) 
equations, fully incorporates the effects of radiation emission, intrachannel in¬ 
teraction, and temporal inhomogeneities. Thus, in order to check whether sta¬ 
ble long-distance multichannel transmission can indeed be realized by a proper 
choice of linear amplifier gain, it is important to carry out numerical simulations 
with the full coupled-NLS model. 

In the current paper, we take on this important task. For this purpose, we 
employ perturbed coupled-NLS models, which take into account the effects of 
second-order dispersion, Kerr nonlinearity, delayed Raman response, and fre¬ 
quency dependent linear gain-loss. We perform numerical simulations with the 
coupled-NLS models with 2 < < 4 frequency channels for two main trans- 
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mission setups. In the first setup, the soliton sequences propagate through a 
single optical fiber, while in the second setup, the sequences propagate through 
a waveguide coupler. We then analyze the simulations results in comparison 
with the predictions of the predator-prey model of Ref. [13], looking for pro¬ 
cesses leading to transmission stabilization and destabilization. The coupled- 
NLS simulations for single-fiber transmission show that at short-to-intermediate 
distances soliton amplitudes exhibit stable oscillatory dynamics, in excellent 
agreement with the predator-prey model’s predictions. These results mean that 
radiation emission and intrachannel interaction effects can indeed be neglected 
at short-to-intermediate distances. However, at larger distances, we observe 
transmission destabilization due to formation of radiative sidebands, which is 
caused by the effects of Kerr nonlinearity on interchannel soliton collisions. We 
also find that the radiative sidebands for the jth soliton sequence form near the 
frequencies /3fc(z) of the solitons in the neighboring frequency channels. Addi¬ 
tionally, we find that the presence of frequency dependent linear gain-loss in 
a single fiber leads to a moderate increase in the distance along which stable 
transmission is observed. The limited enhancement of transmission stability in 
a single fiber is explained by noting that in this case one cannot employ strong 
linear loss at the frequencies of the propagating solitons, and therefore, one 
cannot efficiently suppress the formation of the radiative sidebands. 

A stronger enhancement of transmission stability might be achieved in a 
nonlinear waveguide coupler, consisting of N nearby waveguides. Indeed, in 
this case one might expect to achieve a more efficient suppression of radiative 
sideband generation by employing relatively strong linear loss outside of the 
central amplification frequency interval for each of the N waveguides in the 
waveguides coupler. To test this prediction, we carry out numerical simula¬ 
tions with the coupled-NLS model for propagation in the waveguide coupler. 
The coupled-NLS simulations show that transmission stability and the valid¬ 
ity of the predator-prey model’s predictions in the waveguide coupler system 
are extended to significantly larger distances compared with the distances in 
the single-fiber system. Furthermore, the simulations for the waveguide coupler 
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show that no radiative sidebands form throughout the propagation. Based on 
these observations we conclude that the enhanced transmission stability in the 
waveguide coupler is a result of the efficient suppression of radiative sideband 
generation by the frequency dependent linear gain-loss in this setup. 

We consider optical solitons as an example for the pulses carrying the in¬ 
formation for the following reasons. First, due to the integrability of the un¬ 
perturbed NLS equation and the shape-preserving property of NLS solitons, 
derivation of the predator-prey model for Raman-induced amplitude dynamics 
is done in a rigorous manner [43]. Second, the soliton stability and shape¬ 
preserving property make soliton-based transmission in broadband fiber optics 
links advantageous compared with other transmission methods imiiiiiiT]. 
Third, as mentioned above, the Raman-induced energy exchange in pulse col¬ 
lisions is similar in linear transmission, conventional soliton transmission, and 
dispersion-managed soliton transmission. Thus, even though pulse dynamics in 
these different transmission systems is different, analysis of soliton-based trans¬ 
mission stabilization and destabilization might give a rough idea about the pro¬ 
cesses leading to stabilization and destabilization of the optical pulse sequences 
in other transmission setups. 

The remainder of the paper is organized as follows. In Section]^ we present 
the coupled-NLS model for 7V-channel transmission in a single fiber together 
with the iV-dimensional predator-prey model for Raman-induced dynamics of 
pulse amplitudes. We then review the results of Ref. |23| for stability analysis of 
the equilibrium states of the predator-prey model. In Section we present the 
results of numerical simulations with the coupled-NLS model for single-fiber 
multichannel transmission and analyze these results in comparison with the 
predictions of the predator-prey model. In Section we present the coupled- 
NLS model for pulse propagation in a nonlinear 7V-waveguide coupler. We 
then analyze the results of numerical simulations with this model and compare 
the results with the predator-prey model’s predictions. Our conclusions are 
presented in Section]^ In Appendix A we discuss the method for determining 
the stable propagation distance from the results of the numerical simulations. 
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2. The propagation model for single-fiber transmission and the predator- 
prey model for amplitude dynamics 

We consider propagation of pulses of light in a single-fiber Wchannel trans¬ 
mission link, taking into account second-order dispersion, Kerr nonlinearity, 
delayed Raman response, and frequency-dependent linear loss or gain. The net 
linear gain-loss is the difference between amplifier gain and fiber loss, where we 
assume that the gain is provided by distributed Raman amplification [48l |49] . 

In addition, we assume that the frequency difference A/3 between adjacent chan¬ 
nels is much larger than the spectral width of the pulses, which is the typical 
situation in many soliton-based WDM systems [HEniEiiisaES]. Under these 
assumptions, the propagation is described by the following system of N per¬ 
turbed coupled-NLS equations [S3]: 

N 

+4^(1 - Sjk)\i^k\^'il^j = 

k^l 

N 

- eR^{l - Sjk) , ( 1 ) 

k=l 

where is proportional to the envelope of the electric field of the jth sequence, 

^ < j < N, z is propagation distance, and t is time |SS]. In Eq. 0 
is the Raman coefficient, gioj) is the net frequency dependent linear gain-loss 
function [56] . tp is the Fourier transform of ip with respect to time, stands 
for the inverse Fourier transform, and 6jk is the Kronecker delta function. The 
second term on the left hand side of Eq. Q describes second-order dispersion 
effects, while the third and fourth terms represent intrachannel and interchannel 
interaction due to Kerr nonlinearity. The first term on the right hand side of 
Eq. Q describes the effects of frequency dependent linear gain or loss, the 
second corresponds to Raman-induced intrachannel interaction, while the third 
and fourth terms describe Raman-induced interchannel interaction. 

The form of the net frequency dependent linear gain-loss function g{uj) is 
chosen so that Raman crosstalk and radiation emission effects are suppressed. 
More specifically, g{ui) is equal to a value gj, required to balance Raman-induced 


amplitude shifts, inside a frequency interval of width W centered about the 
initial frequency of the jth-channel solitons /3j(0), and is equal to a negative 
value gL elsewhere. Thus, g(uj) is given by: 



gj if /3j (0) - fT/2 < a; < /3j(0) + for 1 < j < N, 
gL elsewhere. 


where gL < 0- The width W in Eq. ([^ satisfies 1 < IT < A/3, where 
A/3 = /3j+i(0) — /3j(0) for 1 < j < A — 1. Note that the actual values of the 
gj coefficients are determined by the predator-prey model for collision-induced 
amplitude dynamics, such that amplitude shifts due to Raman crosstalk are 
compensated for by the linear gain-loss. The value of gL is determined such 
that instability due to radiation emission is mitigated. In addition, the value 
of W is determined by the following two factors. First, we require W ^ 1, 
such that the effects of the strong linear loss gL on the soliton patterns and 
on the collision-induced amplitude dynamics are relatively small even at large 
distances. Second, we typically require W < A/3, such that instability due to 
radiation emission is effectively mitigated. In practice, we determine the val¬ 
ues oi gL and W by carrying out numerical simulations with the coupled-NLS 
model while looking for the set of values, which yields the longest stable 
propagation distance. Our simulations show that the optimal gL value is around 
0.5, while W should satisfy W > 10. Figure illustrates a typical linear gain- 
loss function g{uj) for a two-channel system with gi = —0.0045, g 2 = 0.0045, 
gL = —0.5, /3i(0) = —7.5, /32(0) = 7.5, and W = 10. These parameter values 
are used in the numerical simulations, whose results are shown in Fig. EJa). 

In the current paper we study soliton-based transmission systems, and there¬ 
fore the optical pulses in the jth frequency channel are fundamental solitons 
of the unperturbed NFS equation idztpj + dftjjj + 2\^|;j\'^'^pj = 0. The en¬ 
velopes of these solitons are given by 'ij^sj(t,z) = gj exp(ixj)sech{xj), where 

= Vj {t — Uj — 2/3^2:), Xj = Oij + Pj{t — yj) + [gj — /3|) z, and the four pa¬ 
rameters gj , j3j , yj , and aj are related to the soliton amplitude, frequency (and 
group velocity), position, and phase, respectively. The assumption of a large fre- 
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Figure 1: (Color online) An example for the frequency-dependent linear gain-loss function 
5 '(aj), described by Eq. in a two-channel system. 


quency (and group velocity) difference between adjacent channels, means that 
|/3j — ,dfe| 1 for 1 < j 1 < fc < N, and j ^ k. As a result of the 

large group velocity difference, the solitons undergo a large number of interse¬ 
quence collisions. The Raman-induced crosstalk during these collisions can lead 
to signihcant amplitude and frequency shifts, which can in turn lead to severe 
transmission degradation. 

In Ref. [43], we showed that the dynamics of soliton amplitudes in an N- 
channel system can be approximately described by an A^-dimensional predator- 
prey model. The derivation of the predator-prey model was based on the follow¬ 
ing simplifying assumptions. (1) The soliton sequences are deterministic in the 
sense that all time slots are occupied and each soliton is located at the center of 
a time slot of width T, where T 3> 1. In addition, the amplitudes are equal for 
all solitons from the same sequence, but are not necessarily equal for solitons 
from different sequences. This setup corresponds, for example, to return-to-zero 
PSK transmission. (2) The sequences are either (a) inhnitely long, or (b) subject 
to periodic temporal boundary conditions. Setup (a) is an approximation for 
long-haul transmission systems, while setup (b) is an approximation for closed 
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fiber-loop experiments. (3) The linear gain-loss coefficients gj in the frequency 
intervals (/3j(0) — W/2 < w < /3j(0) -|- W/2], defined in Eq. ([^, are determined 
by the difference between distributed amplifier gain and fiber loss. In particular, 
for some channels this difference can be slightly positive, resulting in small net 
gain, while for other channels this difference can be slightly negative, resulting 
in small net loss. (4) Since T 3> 1, the solitons in each sequence are temporally 
well-separated. As a result, intrachannel interaction is exponentially small and 
is neglected. (5) The Raman coefficient and the reciprocaf of the frequency spac¬ 
ing satisfy er l/A/3 ^ 1. Consequently, high-order effects due to radiation 
emission are neglected, in accordance with the analysis of the single-collision 
problem 

By assumptions (l)-(5), the propagating soliton sequences are periodic, and 
as a result, the amplitudes of all pulses in a given sequence undergo the same 
dynamics. Taking into account collision-induced amplitude shifts due to delayed 
Raman response, and single-pulse amplitude changes due to linear gain-loss, we 
obtain the following equation for amplitude dynamics of jth-channel solitons 



N 

9j + C '^{k - j)f{\j - k\)gk 


(3) 


L k=l 

where C = ^erA^/T, and 1 < j < N. The coefficients f{\j — k\) on the right 
hand side of Eq. (|^ are determined by the frequency dependence of the Raman 
gain. In particular, for the commonly used triangular approximation for the 
Raman gain curve mm. in which the gain is a piecewise linear function of the 
frequency, /(|j — A:|) = 1 for 1 < j < and 1 < fc < [IH] . 

In WDM systems it is often desired to achieve steady state transmission, 
in which pulse amplitudes in all channels are equal and constant (independent 
of z) [T]. We therefore look for a steady state of the system ^ in the form 
= g > 0 for 1 < j < A, where g is the desired equilibrium amplitude 
value. This yields the following expression for the gj'. 

N 

9i = -C9Y^{k-j)f{\j-k\). (4) 
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Thus, in order to maintain steady-state transmission with equal amplitudes in 
all channels, high-frequency channels should be overamplihed and low-frequency 
channels should be underamplified, compared with central frequency channels. 
Substituting Eq. 0 into Eq. 0, we obtain the following model for amplitude 
dynamics [IH] : 




N 


(5) 




which has the form of a predator-prey model for N species m- 

The steady states of the predator-prey model 0 with nonzero amplitudes in 
all channels are determined by solving the following system of linear equations: 

N 

- j)f{\j - - ??) = 0, l<j<N. (6) 

k^l 

The trivial solution of Eq. i.e., the solution with = rj > 0 for 

1 < fc < TV, corresponds to steady state transmission with equal nonzero ampli¬ 
tudes. Note that the coefficients (fc — j)f{\j — k\) in Eq. 0 are antisymmetric 
with respect to the interchange of j and k. As a result, for WDM systems with 
an odd number of channels, Eq. 0 has infinitely many nontrivial solutions, 
which correspond to steady states of the predator-prey model 0 with unequal 
nonzero amplitudes. This is also true for WDM systems with an even num¬ 
ber of channels, provided that the Raman gain is described by the triangular 
approximation [43]. 

The stability of all the steady states with nonzero amplitudes, rjj = > 

0,1 < j < was established in Ref. [23j; by showing that the function 


N 


VLiv) = 

i=i L 


Vj 


(eq) . (eq) , / 

rij -f Mn I — 


(eq) ' 


Vj 


(7) 


where rj = ( 771 ,..., ..., 7]^), is a Lyapunov function for the predator-prey 

model 0. This stability was found to be independent of the /(|j — fc|) values, 
i.e., of the specific details of the approximation to the Raman gain curve. Fur¬ 
thermore, since dVL/dz = 0 along trajectories of 0 , rather than dVL/dz < 0 , 
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typical dynamics of the amplitudes r]j(z) for input amplitudes that are off the 
steady state value is oscillatory [43]. This behavior also means that the steady 
states with nonzero amplitudes in all channels are nonlinear centers of Eq. 0 


3. Numerical simulations for single-fiber transmission 

The predator-prey model, described in section is based on several sim¬ 
plifying assumptions, whose validity might break down with increasing number 
of channels or at large propagation distances. In particular, the predator-prey 
model neglects radiation emission and modulation instability, intrasequence in¬ 
teraction, and deviations from the assumed periodic form of the soliton se¬ 
quences. These effects can lead to instabilities and pulse-pattern corruption, 
and also to the breakdown of the predator-prey model description (see, for ex¬ 
ample Refs, im [45l [46] , for the case of crosstalk induced by nonlinear gain or 
loss). In contrast, the coupled-NLS model 0 provides a fuller description of the 
propagation, which includes all these effects. Thus, in order to check the predic¬ 
tions of the predator-prey model ([^ for stable dynamics of soliton amplitudes 
and the possibility to realize stable long-distance multichannel soliton-based 
transmission, it is important to carry out numerical simulations with the full 
coupled-NLS model. 

In the current section, we first present numerical simulations with the sys¬ 
tem 0 without the Raman and the linear gain-loss terms. We then present 
a comparison between simulations with the full coupled-NLS model Q with 
the Raman term and the linear gain-loss profile ([^ and the predictions of the 
predator-prey model ([^ for collision-induced amplitude dynamics. We con¬ 
clude the section by analyzing pulse-pattern deterioration at large distances, as 
observed in the full coupled-NLS simulations. 

The coupled-NLS system Q is numerically solved using the split-step method 
with periodic boundary conditions [T] . The use of periodic boundary conditions 
means that the numerical simulations describe pulse dynamics in a closed fiber 
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loop. The initial condition is in the form of N periodic sequences of 2 J solitons 
with initial amplitudes r]j{0), initial frequencies /3j(0), and initial zero phases: 


.7-1 




k=-J 


7?j(0) exp{if3j{0)[t — {k + 1/2)T — (5^-]} 
cosh{r7j(0)[t — (fc + 1/2)T — ’ 


( 8 ) 


where 1 < j < The coefficients 6j in Eq. ([^ correspond to the initial position 
shift of the pulses in the jth sequence relative to pulses located at (fc + l/2)T for 
—J<k<J—l. We simulate multichannel transmission with two, three, and 
four channels and two solitons in each channel. Thus, 2 < < 4 and J = 1 are 

used in our numerical simulations. To maximize the stable propagation distance, 
we choose /3i(0)=-/32(0) for a two-channel system; /3i(0)=-,53(0), /32(0) = 0 
for a three-channel system; and /3i(0)=-/34(0), /32{0)=-(33{0) for a four-channel 
system. In addition, we take 6j = {j — 1)T/N for 1 < j < A^. These choices are 
based on extensive numerical simulations with Eq. 0 and different values of 
Pj{0) and Sj. 

In the numerical simulations, we consider as a concrete example transmission 
at a bit-rate B = 12.5 Gb/s per channel with the following physical parameter 
values inn]. The pulse width and time slot width are r = 5 ps and T = 80 
ps, and the frequency spacing is taken as At/ = 0.48 THz for N = 2,3, and 4 
channels. Thus, the total bandwidth of the system is smaller than 13.2 THz, 
and all channels lie within the main body of the Raman gain curve. The values 
of the dimensionless parameters for this system are er = 0.0012, T = 16, and 
A/3 = 15 for N = 2,3,4. Assuming ^2 = —4 ps^km“^ and 7 = 4 W“^km“^ for 
the second-order dispersion and Kerr nonlinearity coefficients, the soliton peak 
power is Pq = 40 mW. Tables 1 and 2 summarize the values of the dimensionless 
and dimensional physical parameters used in the simulations. In these tables, 
W and W stand for the dimensionless and dimensional width of the linear gain- 
loss function g(w) in Eq. 0, while Zg and Xg correspond to the dimensionless 
and dimensional distance along which stable propagation is observed. 

Note that the Kerr nonlinearity terms appearing in Eq. Q are nonperturba- 
tive. Even though these terms are not expected to affect the shape, amplitude. 
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Table 1: The dimensionless parameters 
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and frequency of a single soliton, propagating in an ultralong optical fiber, the 
situation can be very different for multiple soliton sequences, circulating in a 
fiber loop. In the latter case, Kerr-induced effects might lead to radiation emis¬ 
sion, modulation instability, and eventually to pulse-pattern corruption [H]. It 
is therefore important to first analyze the effects of Kerr nonlinearity alone on 
the propagation. For this purpose, we carry out numerical simulations with the 
following coupled-NLS model, which incorporates second-order dispersion and 
Kerr nonlinearity, but neglects delayed Raman response and linear gain-loss: 

N 

idzi^j + +4^(1 - Sjk)\ipk\‘^'ilJj = 0, (9) 

k^l 

where 1 < j < N. The simulations are carried out for two, three, and four 
frequency channels with the physical parameter values listed in rows 2, 4, and 
7 of Table 1. As an example, we present the results of the simulations for the 
following sets of initial soliton amplitudes: ?7i(0) = 0.9, 772 ( 0 ) = 1.05 for TV = 2; 
771 ( 0 ) = 0.9, 772 ( 0 ) = 0.95, 773 ( 0 ) = 1.1 for N = 3] and 771 ( 0 ) = 0.9, 772 ( 0 ) = 0.95, 
773 ( 0 ) = 1.05, 774 ( 0 ) = 1.15 for = 4. We emphasize, however, that similar 
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Table 2: The dimensional parameters 


jYo 

N 

To (ps) 

r(ps) 

W (THz) 

Xs (km) 


Figures 

1 

2 

5 

80 

0.32 

11875 1 


: a 

>1 


:a)-(b) 

2 

2 

5 

80 

0 

6875 


1 

a) 

-(b) 

3 

2 

5 

80 

0.32 

140000 


a). 


:a)-(b) 

4 

3 

5 

80 

0 

6375 


B 

c) 

-(d) 

5 

3 

5 

80 

0.32 

7750 

4 

b), 

5 

i 

:)-(d) 

6 

3 

5 

80 

0.32 

150625 

8 

b), 

9 

( 

=)-(d) 

7 

4 

5 

80 

0 

4250 


i 

e)-(f) 

8 

4 

5 

80 

0.35 

6250 

4 

c), 

5 



9 

4 

5 

80 

0.35 

45000 

8 

c), 

9 



10 

4 

5 

80 

0.48 

22500 
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results are obtained with other choices of the initial soliton amplitudes. The 
numerical simulations are carried out up to a distance Zg, at which instabil¬ 
ity appears. More specifically, we define Zg as the largest distance at which 


the values of the integrals Ij{z) in Eq. (A.3) in Appendix A are still smaller 
than 0.05 for 1 < j < N. The actual value of Zg depends on the values of the 
physical parameters and in particular on the number of channels N. For the 
coupled-NLS simulations with Eq. and the aforementioned initial amplitude 
values, we find Zg-^ = 550 for N = 2, Zg^ = 510 for N = 3, and Zsg = 340 for 
N = 4. Figurej^shows the pulse patterns Zg)\ and their Fourier transforms 
\tpj(u}^ Zg)\ at the onset of instability, as obtained by the numerical solution of 
Eq. ([^. Also shown are the theoretical predictions for the pulse patterns and 
their Fourier transforms at the onset of instability. Figure shows magnified 
versions of the graphs in Fig. [^for small \^j{t,Zg)\ and \tjjj{u}, Zs)\ values. The 
theoretical prediction for \tjjj{t,Zs)\ is obtained by summation over fundamental 
NTS solitons with amplitudes r]j{0), frequencies /3j(0), and positions yj{zg) + kT 


for —J<fc<J—1, which are measured from the simulations (see Appendix 
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[a| . The theoretical prediction for Zs)\ is obtained by taking the Fourier 

transform of the latter sum. As can be seen from Fig. the soliton patterns 
are almost intact at z = for N = 2,3,4. Additionally, the soliton amplitude 
and frequency values are very close to their initial values. Thus, the solitons 
propagate in a stable manner up to the distance Zs- However, an examination 
of Figs. 3(a), 3(c), and 3(e) reveals that the soliton patterns are in fact slightly 
distorted at Zg, and that the distortion appears as fast oscillations in the soli¬ 
tons tails. Furthermore, as seen in Figs. 3(b), 3(d), and 3(f), the distortion is 
caused by resonant generation of radiaitive sidebands, where the largest side¬ 
bands for the jth soliton sequence form at frequencies /3j_i(0) and/or /3j+i(0) of 
the neighboring soliton sequences. In addition, the amplitudes of the radiative 
sidebands increase as the number of channels increases (see also Ref. [T7] for 
similar behavior), and as a result, the stable propagation distance Zg decreases 
with increasing JV. The growth of radiative sidebands and pulse distortion with 
increasing z eventually leads to the destruction of the soliton sequences. We 
point out that when each soliton sequence propagates through the fiber on its 
own, no radiative sidebands develop and no instability is observed up to dis¬ 
tances as large as z = 20000 m- The latter finding is also in accordance with 
results of single-channel soliton transmission experiments, which demonstrated 
stable soliton propagation over distances as large as 10® km [60j. Based on these 
observations we conclude that transmission instability in the multichannel opti¬ 
cal fiber system is caused by the Kerr-induced interaction in interchannel soliton 
collisions, that is, it is associated with the terms 2|'0fcP'!/’j Eq. Q. 

We now take into account the effects of delayed Raman response and fre¬ 
quency dependent linear gain-loss on the propagation. Our first objective is to 
check the validity of the predator-prey model’s predictions for collision-induced 
dynamics of soiton amplitudes in the presence of delayed Raman response. For 
this purpose, we carry out numerical simulations with the full coupled-NLS 
model Q with the linear gain-loss function ([^ for two, three, and four fre¬ 
quency channels with the physical parameter values listed in rows 1,5, and 8 of 
Table 1. To enable comparison with the results presented in Figs. [^and[^ we 
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Figure 2: (Color online) The pulse patterns at the onset of transmission instability Zs)\ 

and their Fourier transforms Zs)\ for two-channel [(a)-(b)], three-channel [(c)-(d)], and 

four-channel [(e)-(f)] transmission in the absence of delayed Raman response and linear gain- 
loss. The physical parameter values are listed in rows 2, 4, and 7 of Table 1. The stable 
transmission distances are Zsi = 550 for N = 2, Zs 2 = 510 for N = 3, and Zs^ = 340 for 
= 4. The solid-crossed red curve [solid red curve in (a)], dashed green curve, solid blue 
curve, and dashed-dotted magenta curve represent Zs)\ with j = 1,2, 3,4, obtained by 

numerical simulations with Eq. §. The red circles, green squares, blue up-pointing triangles, 
and magenta down-pointing triangles represent |'0j(t,2s)| with j = 1,2,3,4, obtained by the 
simulations. The brown diamonds, gray left-pointing triangles, black right-pointing triangles, 
and orange stars represent the theoretical prediction for Zs)\ or Zs)\ with j = 

1, 2, 3,4, respectively. 
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Figure 3: (Color online) Magnified versions of the graphs in Fig. [^for small Zs)\ and 

Zs)\ values. The symbols are the same as in Fig. 
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discuss the results of simulations with the same sets of initial soliton amplitudes 
as the ones used in Figs. [^and|^ The numerical simulations are carried out up 
to the onset of transmission instability, which occurs at = 950 for N = 2, 
at Zs^ = 620 for N = 3, and at Zg^ = 500 for = 4. The z dependence of soli¬ 
ton amplitudes obtained by numerical solution of the full coupled-NLS model 
Q with the gain-loss ([^ is shown in Fig. along with the prediction of the 
predator-prey model <§■ In all three cases the soliton amplitudes oscillate about 
their equilibrium value 77 = 1, i.e., the dynamics of soliton amplitudes is stable 
up to the distance Zg. Furthermore, the agreement between the coupled-NLS 
simulations and the predator-prey model’s predictions is excellent throughout 
the propagation. Thus, our coupled-NLS simulations validate the predictions 
of the predator-prey model ([^ for collision-induced amplitude dynamics in the 
presence of delayed Raman response at distances Q < z < Zg. This is a very 
important observation, because of the major simplifying assumptions that were 
made in the derivation of the model ([^. In particular, we conclude that the ef¬ 
fects of radiation emission, modulation instability, intrachannel interaction, and 
other high-order perturbations can indeed be neglected for distances smaller 
than Zg. 

Further insight about transmission stability and about the processes leading 
to transmission destabilization is gained by an analysis of the soliton patterns 
at the onset of instability. Figure shows the pulse patterns at the onset of 
instability \il}j{t^Zs)\ and their Fourier transforms Zg)\., obtained by the 

numerical simulations that are described in the preceding paragraph. The the¬ 
oretical predictions for the pulse patterns and their Fourier transforms, which 
are calculated in the same manner as in Fig. are also shown. In addition. 
Fig. § shows magnified versions of the graphs in Fig. [^for small \'iljj{t, Zg)\ and 
\^pj(u}., Zs)| values. We observe that the soliton patterns are almost intact at Zg. 
Based on this observation and the observation that dynamics of soliton ampli¬ 
tudes is stable for 0 < z < we conclude that the multichannel soliton-based 
transmission is stable at distances smaller than Zg. However, as seen in Figs. 
[^a),gc), and[^ e), the soliton patterns are actually slightly distorted at Zg, 
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Figure 4: (Color online) The 2 : dependence of soliton amplitudes rjj for two-channel (a), three- 
channel (b), and four-channel (c) transmission in the presence of delayed Raman response 
and the frequency dependent linear gain-loss §. The physical parameter values are listed 
in rows 1, 5, and 8 of Table 1. The red circles, green squares, blue up-pointing triangles, 
and magenta down-pointing triangles represent 771 ( 2 ), 772 ( 2 ), 773 ( 2 ), and 774 ( 2 ), obtained by 
numerical solution of Eqs. 0 and §. The solid brown, dashed gray, dashed-dotted black, 
and solid-starred orange curves correspond to 771 ( 2 ), 772 ( 2 ), 773 ( 2 ), and 774 ( 2 ), obtained by the 
predator-prey model 0 - 
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and the distortion appears as fast oscillations in the solitons tails. Moreover, 
as seen in Figs. [^b),|^d), and|^f), pulse-pattern distortion is caused by res¬ 
onant formation of radiative sidebands, where the largest sidebands for the jth 
sequence form near the frequencies j3j-i{z) and/or j5j+i[z) of the neighboring 
soliton sequences. Thus, the mechanisms leading to deterioration of the soli- 
ton sequences in multichannel transmission in a single fiber in the presence of 
delayed Raman response and linear gain-loss are very similar to the ones ob¬ 
served in the absence of delayed Raman response and linear gain-loss. This also 
indicates that the dominant cause for transmission destabilization in the full 
coupled-NLS simulations for multichannel transmission in a single fiber is due 
to the effects of Kerr-induced interaction in interchannel soliton collisions. As 
explained earlier, the latter effects are represented by the 2|'0fcp'(/;j terms in Eq. 

0 . 

It is interesting to note that the stable propagation distances Zs in the pres¬ 
ence of delayed Raman response and linear gain-loss are larger compared with 
the distances obtained in the absence of these two processes by factors of 1.7 for 
iV = 2, 1.2 for iV = 3, and 1.5 for N = A. We attribute this moderate increase 
in Zs values to the introduction of frequency dependent linear gain-loss with 
strong loss ql outside the frequency intervals /3j(0) — W/2 < uj < /3j(0) -I- lE/2, 
where 1 < j < A^, which leads to partial suppression of radiative sideband 
generation. However, the suppression of radiative instability in a single fiber is 
quite limited, since the radiative sidebands for a given sequence form near the 
frequencies f3k{z) of the other soliton sequences. As a result, in a single fiber, 
one cannot employ strong loss at the latter frequencies, as this would lead to 
the decay of the propagating solitons. Better suppression of radiative instability 
and significantly larger Zg values can be realized in nonlinear waveguide couplers 
with frequency dependent linear gain-loss. This subject is discussed in detail in 
Section |4l 

We now turn to discuss the later stages of pulse pattern deterioration, i.e., 
the evolution of the soliton sequences in a single fiber for distances z > Zs- As 
a concrete example, we discuss the four-channel setup considered in Figs. Qc), 
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Figure 5: (Color online) The pulse patterns at the onset of transmission instability \'ipj[t^Zs)\ 
and their Fourier transforms Zs)\ for two-channel [(a)-(b)], three-channel [(c)-(d)], and 

four-channel [(e)-(f)] transmission in the presence of delayed Raman response and the linear 
gain-loss The physical parameter values are listed in rows 1, 5, and 8 of Table 1. The stable 
transmission distances are Zs^^ = 950 for N = 2, Zs^ = 620 for = 3, and Zsq = 500 for = 4. 
The solid-crossed red curve [solid red curve in (a)], dashed green curve, solid blue curve, and 
dashed-dotted magenta curve represent Zs)\ with j = 1,2,3, 4, obtained by numerical 

simulations with Eqs. 0 and 0. The red circles, green squares, blue up-pointing triangles, 
and magenta down-pointing triangles represent Zs)\ with j = 1,2,3,4, obtained by the 

simulations. The brown diamonds, gray left-pointing triangles, black right-pointing triangles, 
and orange stars represent the theoretical prediction for \'ipj{t^Zs)\ or \'ipj(u}, Zs)\ with j = 
1, 2, 3,4, respectively. 
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Figure 6: (Color online) Magnified versions of the graphs in Fig. [^for small Zs)\ and 

2 : 5)1 values. The symbols are the same as in Fig. 
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[^e), and[^f), for which Zgg = 500 [SH]- Figures [^a) and[^b) show the pulse 
patterns \ipj(t,z)\ and their Fourier transforms \ipj(uj,z)\ at z = 600 > Zsg, 
as obtained by numerical solution of Eqs. 0 and 1^. It is seen that the 
largest sidebands at z = 600 form near the frequency Psiz) for the j = 2 
and j = 4 soliton sequences, and near the frequencies /32{z) and /34(z) for the 
j = 3 soliton sequence. These larger sidebands lead to significantly stronger 
pulse distortion at z = 600 compared with Zgj, = 500. In particular, at z = 
600, the j = 3 pulse sequence is strongly distorted, where the distortion is in 
the form of fast oscillations in the main body of the solitons. In contrast, at 
Zgj, = 500, the j = 3 sequence is only weakly distorted, and the distortion is 
in the form of fast oscillations, which are significant only in the solitons tails. 
Additionally, radiation emitted by the solitons in the j = 2, j = 3, and j = 4 
frequency channels develops into small pulses at z = 600. The largest radiation- 
induced pulses are generated due to radiation emitted by solitons in the j = 4 
channel near the frequency Psiz). Figures [^c) and 0d) show a comparison 
of the shape and Fourier transform of the latter pulses with the shape and 
Fourier transform expected for a single NTS soliton with the same amplitude 
and frequency. It is clear that these radiation-induced pulses do not posses the 
soliton form. Similar conclusion holds for the other radiation-induced pulses. 
The amplitudes of the radiative sidebands generated by the j = 2, j = 3, and 
J = 4 pulse sequences continue to increase with increasing propagation distance 
and this leads to further pulse pattern degradation. Indeed, as seen in Fig. 
0f)> at z = 650, the radiative sidebands generated by the j = 4 sequence near 
/33 (z) and by the j = 3 sequence near /3i{z) are comparable in magnitude to 
the Fourier transforms of the j = 3 and j = 4 pulse sequences, respectively. 
Additional strong radiative sidebands are observed for the j = 2 sequence near 
frequencies Psiz) and P 4 {z) and for the j = 4 sequence near frequency P 2 iz). As 
a result, the j = 2, j = 3, and j = 4 pulse sequences are strongly degraded due to 
pulse distortion at z = 650. More specifically, distortion due to fast oscillations 
in both the main body and the tail of the pulses is observed for these three 
pulse sequences [see Fig. [^e)]. In addition, the number and amplitudes of 
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the radiation-induced pulses are much larger at z = 650 compared with the 
corresponding number and amplitudes of these pulses at z = 600. 

4. Nonlinear waveguide coupler transmission 

The results of the numerical simulations in Section |3] show that in a sin¬ 
gle fiber, radiative instabilities can be partially mitigated by employing the 
frequency dependent linear gain-loss ([^. However, as described in Section]^ 
suppression of radiation emission in a single hber is still quite limited, and gen¬ 
eration of radiative sidebands leads to severe pulse pattern degradation at large 
distances. The limitation of the single-hber setup is explained by noting that 
the radiative sidebands for each pulse sequence form near the frequencies /3fc(z) 
of the other pulse sequences. As a result, in a single fiber, one cannot employ 
strong loss at or near the frequencies f3k{z), as this would lead to the decay of 
the propagating pulses. It is therefore interesting to look for other waveguide 
setups that can significantly enhance transmission stability. A very promising 
approach for enhancing transmission stability is based on employing a nonlinear 
waveguide coupler, consisting of N very close waveguides In this case each 
pulse sequence propagates through its own waveguide and each waveguide is 
characterized by its own frequency dependent linear gain-loss function gj{uj, z) 
m- This enables better suppression of radiation emission, since the linear gain- 
loss of each waveguide can be set equal to the required gj value within a certain 
z-dependent bandwidth {j3j{z) — W/2, l3j{z) + W/2] around the central frequency 
Pj (z) of the solitons in that waveguide, and equal to a relatively large negative 
value gL outside of that bandwidth. This leads to enhancement of transmission 
stability compared with the single fiber setup, since generation of all radiative 
sidebands outside of the interval i/3j(z) — Wj2,/3j(z) + IT/2] is suppressed by 
the relatively strong linear loss gL- 

In the current section, we investigate the possibility to significantly enhance 
transmission stability in multichannel soliton-based systems by employing N- 
waveguide couplers with frequency dependent linear gain-loss. The enhanced 
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Figure 7: (Color online) The pulse patterns \'il)j{t^z)\ and their Fourier transforms \ipj{LO^z)\ 
qX z = 600 [(a) and (b)] and at 2 = 650 [(e) and (f)], obtained by numerical simulations with 
Eqs. 0 and § for the four-channel system considered in Figs. [^and[^ The symbols in 
(a), (b), (e), and (f) are the same as in Fig. (c) The shape of a radiation-induced pulse 
2 ) 1 , generated due to radiation emitted by the j = A soliton sequence at 2 = 600. The 
magenta down-pointing triangles represent the numerically obtained (i; ^)\ with 2 = 600, 

while the solid orange curve corresponds to the shape of a single NLS soliton with the same 
amplitude, (d) The Fourier transform of the radiation-induced pulse in (c). The magenta 
down-pointing triangles represent the numerically obtained (oj, 2 )| with 2 = 600, while 

the solid orange curve corresponds to the Fourier transform of a single NLS soliton with the 
same amplitude and central frequency. 
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transmission stability is also expected to enable observation of the stable os¬ 
cillatory dynamics of soliton amplitudes, predicted by the predator-prey model 
along significantly larger distances compared with the distances observed in 
single-fiber transmission. Similar to the single-fiber setup considered in Section 
we take into account the effects of second-order dispersion, Kerr nonlinearity, 
delyaed Raman response, and linear gain-loss. The main difference between 
the waveguide coupler setup and the single-fiber setup is that the single linear 
gain-loss function g{uj) of Eq. ([^ is now replaced by N 2 -dependent linear 
gain-loss functions gj{uj,z), where 1 < j < fV. Thus, the propagation of the 
pulse sequences through the waveguide coupler is described by the following 
coupled-NLS model: 

N 

idz'ipj + + 4^(1 - 6jk)\'ipkfipj = 

k^l 

N 

- Sjk) , ( 10 ) 

k^l 


where 1 < j < The linear gain-loss function of the jth waveguide gj{uj,z), 


(11) 


appearing on the right-hand side of Eq. (10), is defined by: 

j g, if P, (z) -Wl2<iv< P, (z) + lE/2, 

gj[uj,z) = 

1 if w < Pj{z) - lE/2, or w > Pj{z) + lE/2, 

where the gj coefficients are determined by Eq. (j^, the 2 dependence of the 
frequencies Pj{z) is determined from the numerical solution of the coupled- 


NLS model (10), and gL < 0. Notice the following important properties of 


the gain-loss ( |11[ ). Eirst, the gain-loss gj inside the central frequency interval 
{Pj{z) — W/2,Pj{z) + lE/2] is expected to compensate for amplitude shifts due 
to Raman crosstalk and by this, lead to stable oscillatory dynamics of soliton 
amplitudes. Second, the relatively strong linear loss g^ outside the interval 
{Pj{z) — W/2, Pj{z) + lE/2] should enable efficient suppression of radiative side¬ 
band generation for any frequency outside of this interval. Third, the end points 
of the central frequency interval are shifting with 2 , such that the interval is 
centered around Pj{z) throughout the propagation. This shifting of the central 
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amplification interval is introduced to compensate for the significant Raman- 
induced frequency shifts experienced by the solitons during the propagation 
[52] • The combination of the three properties of gj{uj, z) should lead to a signifi¬ 
cant increase of the stable transmission distances in the nonlinear iV-waveguide 
coupler compared with the single-fiber system considered in Section]^ As a re¬ 
sult, one can expect that the stable oscillatory dynamics of soliton amplitudes, 
predicted by the predator-prey model ([^ , will also hold along significantly larger 
distances. 

In order to check whether the A^-waveguide coupler setup leads to enhance¬ 


ment of transmission stability, we numerically solve Eq. (10) with the gain-loss 


(11) for two, three, and four channels. The comparison with results obtained 
for single-fiber transmission is enabled by using the same values of the phys¬ 
ical parameters that were used in Figs. Ba),|4l:b) and Be) . The numerical 
simulations are carried out up to the onset of transmission instability, which 
occurs at Zg., = 11200 for N = 2, Zgg = 12050 for iV = 3, and Zgg = 3600 for 

= 4. Figure shows the z dependence of soliton amplitudes as obtained 
by the coupled-NLS simulations along with the prediction of the predator-prey 
model (B. It is seen that the amplitudes exhibit stable oscillations about the 
equilibrium value rj = 1 for N = 2,3, and 4. Furthermore, the agreement be¬ 
tween the coupled-NLS simulations and the predictions of the predator-prey 
model are excellent throughout the propagation. Thus, both transmission sta¬ 
bility and the validity of the predator-prey model’s predictions are extended to 
distances that are larger by factors of 11.8 for N = 2, 19.4 for N = 3, and 7.2 
for A^ = 4 compared with the distances obtained with the single-fiber WDM 
system in Section]^ 

Further insight into the enhanced transmission stability in waveguide cou¬ 
plers is gained by analyzing the pulse patterns at the onset of instability I'lpjit, Zs)| 
and their Fourier transforms Zs)\. Figure|^shows the results obtained by 

numerical solution of Eqs. (10) and 0 together with the theoretical predic¬ 
tion. Figure shows magnified versions of the graphs in Fig. for small 
\ipj{t,Zs)\ and \ipj{uj,Zs)\ values. It is seen that the soliton patterns are almost 
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1.2 


(a) 



Figure 8: (Color online) The z dependence of soliton amplitudes rjj in a two-channel (a), a 
three-channel (b), and a four-channel (c) nonlinear waveguide coupler with linear gain-loss 
The values of the physical parameters are the same as the ones used in Fig. The 
red circles, green squares, blue up-pointing triangles, and magenta down-pointing triangles 
represent r)i{z), r]2(z), ^73(2), and 774(2:) obtained by numerical solution of the coupled-NLS 
model ( |10[ l with the gain-loss {nj. The solid brown, dashed gray, dashed-dotted black, 
and solid-starred orange curves correspond to 771(2:), 772(2), 773(2), and 774(2) obtained by 
the predator-prey model 
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intact at z = Zg, in accordance with our conclusion about transmission stability 


for 0 < z < Zg. Moreover, as seen from Fig. 10 no radiative sidebands and no 
fast oscillations in the solitons shapes are observed at z = Zg. Instead, soliton 
distortion at z = Zg is due to slow variations in the shape at the pulse tails, 
formation of small radiative pulses, and position shifts of the solitons from the 
same sequence relative to one another. Thus, instability due to formation of 
radiative sidebands is completely suppressed in A^-waveguide couplers with the 
gain-loss 0- This finding explains the significant increase in the values of 
the stable propagation distance in A^-waveguide coupler transmission compared 
with the distances achieved in single-fiber transmission in Section 

We note that the frequency shifts experienced by the propagating solitons 
at large propagation distances are quite large for both the single-fiber sys¬ 
tems of Section and the A^-waveguide coupler systems of the current sec¬ 
tion. For example, the total frequency shifts measured at Zgg = 3600 from 
the coupled-NLS simulations for the four-channel waveguide coupler of Fig. 
§c) are A/3i(zgJ = -7.354, ^h{zs,) = -7.296, = -7.284, and 

A/34(zgg) = —7.2204. These Raman-induced frequency shifts make the gain-loss 
functions with fixed frequency intervals [such as the function in Eq. ([^] less ef¬ 
fective in stabilizing soliton amplitude dynamics. In order to compensate for the 
effects of these frequency shifts, a shifting of the central amplification interval 
was introduced into the gain-loss function 0- We now complete the analysis 
of transmission stabilization in the waveguide coupler, by evaluating the im¬ 
pact of shifting of the amplification interval in the gain-loss function 0. For 
this purpose, we consider the following alternative gain-loss functions with fixed 
amplification intervals: 


= 


if 13j (0) - IF/2 <oj<Pj (0) + IF/2, 


( 12 ) 


9L if w < ^i(O) - IF/2, or w > /3j (0) + IF/2, 
where 1 < j < N. We carry out numerical simulations with Eq. (|I0[) and the 


gain-loss functions (121 with IF = 15 for a four-channel waveguide coupler. The 


values of the other physical parameters are the same as the ones used in Fig. 
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Figure 9: (Color online) The pulse patterns at the onset of instability 2 :s)| and their 

Fourier transforms Zs)\ for the two-channel [(a)-(b)], the three-channel [(c)-(d)], and 

the four-channel [(e)-(f)] waveguide couplers of Fig. The final distances are Zsj = 11200 
in (a)-(b), Zs^ = 12050 in (c)-(d), and Zsg = 3600 in (e)-(f). The solid-crossed red curve 
[solid red curve in (a)], dashed green curve, solid blue curve, and dashed-dotted magenta 
curve represent Zs)\ with j = 1,2,3, 4, obtained by numerical simulations with Eqs. 

and The red circles, green squares, blue up-pointing triangles, and magenta down¬ 
pointing triangles represent Zs)\ with j = 1,2, 3, 4, obtained by the simulations. The 

brown diamonds, gray left-pointing triangles, black right-pointing triangles, and orange stars 
represent the theoretical prediction for \il)j(t,Zs)\ or Zs)\ with j = 1,2, 3,4, respectively. 
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Figure 10: (Color online) Magnified versions of the graphs in Fig. [^for small \ipj{t^Zs)\ and 
Zs)\ values. The symbols are the same as in Fig. 
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[^c). The simulations are carried out up to the onset of transmisison instability 
at Zsjg = 1800. Figure[TT]shows the z dependence of soliton amplitudes obtained 
by these coupled-NLS simulations along with the prediction of the predator-prey 
model ([^. We observe stable oscillatory dynamics and good agreement with the 
predator-prey model’s prediction throughout the propagation. We note that the 
stable propagation distance Zs^, = 1800 is larger by a factor of 3.6 compared with 
the stable propagation distance for the corresponding single-fiber system [see 
Fig. s c)], but smaller by a factor of 0.5 compared with the distance obtained 
with the waveguide coupler and the gain-loss Based on this comparison we 
conclude that the introduction of shifting of the central amplification interval 
does lead to an enhancement of transmission stability. On the other hand, we 
also observe that even in the absence of shifting of the amplihcation interval, 
the waveguide coupler setup enables stable propagation along significantly larger 
distances compared with the single-fiber systems considered in Section 

We conclude this section by summarizing the dependence of the stable prop¬ 
agation distance Zs on the number of channels N in the main transmission setups 
considered in the paper. Figure shows the Zg values obtained by numerical 
simulations for single-hber transmission with and without the effects of delayed 
Raman response and the linear gain-loss ([^. The Zg values obtained by the 
simulations for Wwaveguide coupler transmission with the linear gain-loss (111 
are also shown. We observe that in single-hber transmission, the introduction 
of the gain-loss ([^ leads to a moderate increase in the Zg values, despite of the 
presence of delayed Raman response. Moreover, the Zg values obtained in N- 
waveguide coupler transmission are signihcantly larger than the ones obtained in 
single-hber transmission. As explained earlier, the enhanced transmission sta¬ 
bility in waveguide couplers can be attributed to the more efficient suppression 
of radiative sideband generation by the linear gain-loss 
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Figure 11 : (Color online) The 2: dependence of soliton amplitudes rjj in a four-channel nonlin¬ 
ear waveguide coupler with linear gain-loss © and W = 15 . The values of the other physical 
parameters are the same as the ones used in Fig. [^(c). The red circles, green squares, blue 
up-pointing triangles, and magenta down-pointing triangles represent 771(2), 772(2), 773(2), and 
774(2) obtained by numerical solution of the coupled-NLS model with the gain-loss 
The solid brown, dashed gray, dashed-dotted black, and solid-starred orange curves correspond 
to 771(2), 772(2), 773(2), and 774(2) obtained with the predator-prey model 0. 
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Figure 12 : (Color online) The dependence of the stable propagation distance Zs on the number 
of channels N in different transmission setups. The red circles represent the values obtained 
for single-fiber transmission in the presence of delayed Raman response and the linear gain- 
loss by numerical solution of Eqs. 0 and 0. The purple triangles represent the values 
obtained for single-fiber transmission in the absence of delayed Raman response and linear 
gain-loss by numerical solution of Eq. The green squares represent the values obtained 

for W-waveguide coupler transmission by numerical solution of Eqs. and E 3 - 
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5. Conclusions 


We investigated transmission stabilization and destabilization and dynam¬ 
ics of pulse amplitudes induced by Raman crosstalk in multichannel soliton- 
based optical waveguide systems with N frequency channels. We considered 
two main transmission setups. In the first setup, the N soliton sequences prop¬ 
agate through a single optical fiber, while in the second setup, the sequences 
propagate through a waveguide coupler, consisting of N close waveguides. We 
studied the transmission by performing numerical simulations with coupled- 
NLS models, which take into account second-order dispersion, Kerr nonlinear¬ 
ity, delayed Raman response, and frequency dependent linear gain-loss. The 
simulations were carried out for two, three, and four frequency channels in both 
single-fiber and waveguide coupler setups. The results of the coupled-NLS sim¬ 
ulations were compared with the predictions of a simplified predator-prey model 
for dynamics of pulse amplitudes [43], which incorporates amplitude shifts due 
to linear gain-loss and delayed Raman response, but neglects radiation emission 
and intrachannel interaction. The predator-prey model in Ref. |43] predicts 
that stable dynamics of soliton amplitudes can be realized by a suitable choice 
of amplifier gain in different frequency channels. One major goal of our study 
was to validate this prediction. A second major goal was to characterize the 
processes that lead to transmission destabilization and to develop waveguide 
setups, which lead to significant enhancement of transmission stability. 

We first studied soliton-based multichannel transmission in a single fiber in 
the absence of delayed Raman response and linear gain-loss. We found that in 
this case, transmission destabilization is caused by resonant formation of radia¬ 
tive sidebands, where the largest sidebands for the jth soliton sequence form at 
frequencies /3j_i(0) and / or /3j_|_i(0) of the solitons in the neighboring frequency 
channels. Additionally, the amplitudes of the radiative sidebands increase with 
increasing number of channels N, and as a result, the stable propagation dis¬ 
tances decrease with increasing N. Furthermore, the stable propagation dis¬ 
tances obtained in our numerical simulations are significantly smaller compared 
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with the distances achieved in Ref. m for single-channel transmission with the 
same values of the physical parameters. Based on these findings we conclude 
that destabilization of multichannel soliton-based transmission in a single fiber 
is caused by Kerr-induced interaction in interchannel soliton collisions. 

We then carried out numerical simulations for multichannel transmission in 
a single fiber, taking into account the effects of delayed Raman response and 
frequency dependent linear gain-loss. We assumed that the gain-loss function 
g{io) for single-fiber transmission is given by Eq. ([^. That is, g(w) is equal 
to the constants gj, determined by the predator-prey model, in frequency in¬ 
tervals (/3j(0) — W/2,/3j(0) -I- W/2] of constant width W centered about the 
initial soliton frequencies PjiO), and is equal to a negative value gL outside of 
these intervals. Numerical simulations with the full coupled-NLS model showed 
that at distances smaller than the stable propagation distance Zg, soliton am¬ 
plitudes exhibit stable oscillatory dynamics, in excellent agreement with the 
predictions of the predator-prey model of Ref. [33]. These findings are very 
important because of the major simplifying assumptions made in the derivation 
of the predator-prey model. In particular, based on these findings, we conclude 
that the effects of radiation emission and intrachannel interaction can indeed 
be neglected at distances smaller than Zg. However, at distances z cz Zg, we 
observed transmission destabilization due to formation of radiative sidebands. 
The destabilization process is very similar to destabilization in the absence of 
delayed Raman response and linear gain-loss, i.e., the largest sidebands for the 
jth soliton sequence form at frequencies Pj-i{z) and / or /3j+i(z) of the solitons 
in the neighboring frequency channels. At distances larger than Zg, the contin¬ 
ued growth of the radiative sidebands leads to fast oscillations in the main body 
of the solitons and to generation of new pulses, which do not possess the soliton 
sech form. As a result, the pulse patterns at this late stage of the propagation 
are strongly distorted. 

We note that the stable propagation distances Zg for single-fiber multichan¬ 
nel transmission in the presence of delayed Raman response and linear gain-loss 
are larger compared with the distances obtained in the absence of these pro- 
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cesses. We attribute this increase in Zs values to the introduction of frequency 
dependent linear gain-loss with relatively strong loss gL outside the frequency 
intervals (/3j(0) — W/2,l3j{0) + W/2], which leads to partial suppression of ra¬ 
diative sideband generation. However, the suppression of radiative instability in 
single-fiber transmission is quite limited, since the radiative sidebands for each 
sequence form near the frequencies Pkiz) of the other soliton sequences. As a 
result, in a single fiber, one cannot employ strong loss at the latter frequencies, 
as this would lead to the decay of the propagating solitons. 

A more promising approach for achieving significant enhancement of trans¬ 
mission stability is based on employing a nonlinear waveguide coupler, con¬ 
sisting of N close waveguides. In this case, each soliton sequence propagates 
through its own waveguide, and each waveguide is characterized by its own fre¬ 
quency dependent linear gain-loss function. We assumed that the linear gain-loss 
for the jth waveguide gj{uj,z) is equal to the constant gj, determined by the 
predator-prey model, inside a z-dependent frequency interval centered about 
the soliton frequency Pjiz), and is equal to a negative value gL outside of this 
interval. This waveguide coupler setup is expected to lead to enhanced transmis¬ 
sion stability, since generation of all radiative sidebands outside of the central 
amplification interval is suppressed by the relatively strong linear loss gL for 
each of the N waveguides. To test this prediction, we carried out numerical 
simulations with a new coupled-NLS model, which takes into account the ef¬ 
fects of second-order dispersion, Kerr nonlinearity, delayed Raman response, 
and the N frequency-dependent linear gain-loss functions gj(uj,z). The sim¬ 
ulations with the new coupled-NLS model showed that transmission stability 
and the validity of the predator-prey model’s predictions in waveguide coupler 
transmission are extended to distances that are larger by factors of 11.8 for two 
channels, 19.4 for three channels, and 7.2 for four channels, compared with the 
distances in the single-fiber system. Additionally, the simulations showed that 
the solitons retain their shape at distances smaller than Zs and no radiative 
sidebands appear throughout the propagation. Based on these observations we 
conclude that transmission stability in the waveguide coupler system is indeed 
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significantly enhanced compared with the single-fiber system. Furthermore, the 
enhanced transmission stability in the waveguide coupler is enabled by the ef¬ 
ficient suppression of radiative sideband generation due to the presence of the 
linear gain-loss gj{uj,z). 

To complete the analysis of transmission stabilization in the waveguide cou¬ 
pler, we evaluated the impact of the shifting of the central amplification in¬ 
tervals of the gain-loss functions gj{LO,z). For this purpose, we replaced each 
gj{uj,z) by a z-independent gain-loss function gj{uj) that is equal to gj inside a 
fixed frequency interval centered about the initial soliton frequency and 

is equal to a negative value g^ outside of this interval. Numerical simulations 
with the coupled-NLS model for a four-channel system showed that the stable 
propagation distance and the distance along which the predator-prey model’s 
predictions are valid are larger by a factor of 3.6 compared with the distance 
achieved in the single-fiber system, but smaller by a factor of 0.5 compared with 
the distance obtained in the waveguide coupler with shifting of the central am¬ 
plification intervals of the linear gain-loss. Based on these findings we conclude 
that the introduction of shifting of the central amplification intervals does lead 
to enhanced transmission stability. On the other hand, even in the absence of 
shifting of the central amplification intervals, the waveguide coupler setup en¬ 
ables stable propagation along significantly larger distances compared with the 
single-fiber setup. 

Appendix A. The method for determining the stable propagation 
distance 

In this Appendix, we present the method that we used for determining the 
value of the stable propagation distance Zg from the results of the coupled- 
NLS simulations. In addition, we present the theoretical predictions for the 
soliton patterns and their Fourier transforms, which were used in the analysis 
of transmission stability. 

We consider propagation of the soliton sequence in the jth frequency channel 
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through an optical waveguide, where the envelope of the electric field of this 
sequence at 2 = 0 is given by Eq. We are interested in the envelope of 
the pulse sequence at distance z. We assume that the pulse sequence is only 
weakly distorted at this distance. Thus, by the standard adiabatic perturbation 
technique for the NLS soliton, one can write the envelope of the electric field of 
the jth sequence at distance z as: z) = iptjit, z) + Vrj(t, z), where iptjit, z) 


is the soliton part and Vrj{t,z) is the radiation part 


Since T ^ 1, 


we neglect intrasequence interaction, in accordance with the assumptions of the 
predator-prey model of Ref. [43] (see also Section]^. Under this assumption, 
the soliton part of the electric field can be expressed as: 


■MU) = ,,(.)e».V. g (A.1) 

^ cosh{r]j[z)[t-yj(z)-kT\} 


k=-J 


where rjj{z) is the amplitude, j5j{z) is the frequency, Oj{z) is the common overall 
phase, yj{z) = Ayj{z) + T/2 + Sj, and Ayj{z) is the common overall position 
shift. The Fourier transform of iptjit, z) with respect to time is given by: 

^tj{u:,z)= (^) sech ~ ^ (A.2) 

( zrij{z) j k=-j 

Our theoretical prediction for the jth pulse pattern at distance z, z)\, is 


calculated by using Eq. (A.l) with values of rjjiz), Pj{z), and yj{z), which are 
measured from the numerical simulations. Similarly, our theoretical prediction 
for the Fourier transform of the jth pulse pattern, z)|, is obtained by 


using Eq. (A.2) with the numerically obtained values of r]j{z) and /3j(z) [65] . 

The method for determining the stable propagation distance is based on a 
comparison of the theoretical predictions for the pulse patterns z)\ with 

the results of the numerical simulation (t, z)| for 1 < j < fV. More 

specifically, we calculate the following normalized integrals, which measure the 
deviation of the numerically obtained pulse patterns from the corresponding 
theoretical predictions: 


i,{z) = iff\z)/i,{z), 


(A.3) 
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where 


if^\z) = 


JT 


~JT 


1/2 




I (num) 
^3 


(t.z) 


dt 


(A.4) 


■ JT 


■ 

j 


2 

dt 

-JT 




(A.5) 


and 1 < j < N. We then define the stable propagation distance Zg as the 
largest distance at which the values of Ij{z) are still smaller than a constant 
C for 1 < j < N. In practice, we used the value C = 0.05 in the numerical 
simulations. We emphasize, however, that the values of the stable propagation 
distance obtained by this method are not very sensitive to the choice of the 
constant C. That is, we found that small changes in the value of C lead to 
small changes in the measured Zg values. 
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